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SUMMARY 



A method for determining lateral stability derivatives from flight 
measurements is obtained by arranging the lateral equations of motion in 
such form that information from each of the three modes of lateral motion 
may be utilized. This method permits determination of all the important 
» derivatives without requiring an estimation of any of the derivatives. 

The results of an error analysis are given to show the effects of errors 
in the measured quantities on the accuracy of each stability derivative 
* for three representative airplanes. 


INTRODUCTION 


Calculated and measured wind-tunnel values of the aerodynamic lateral 
stability derivatives are necessary in the design of any modern airplane. 
Because of the uncertainty of determination of many of the derivatives 
from theory or from wind-tunnel tests, however, flight measurements of 
the derivatives are desirable to check the values assumed in the design 
and to provide a basis for further improvement of the aerodynamic 
characteristics . 

Several methods have been proposed for determining lateral stability 
derivatives from flight tests by analyzing transient or frequency-response 
data. (For example, see refs. 1 to 3*) These methods allow determination 
of all the important derivatives. One method, known as the vector method 
(ref. 3 ), allows some insigjit into the effect of errors in the flight 
measurements on the accuracy of the derivatives. In this method, which 
utilizes data from the Dutch roll mode alone, two of the derivatives must 
be estimated or assumed in order to evaluate the others. 
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The purpose of this paper i s to present a method in which the char- 
acteristics of the spiral and roll-subsidence modes as well as those of 
the Dutch roll mode are utilized in an effort to obtain all the important 
derivatives. In general, separate flight tests are needed to measure 'as 
accurately as possible the characteristics of the three modes. An error 
analysis is made to show the accuracy of flight measurements required to 
produce a desired accuracy of each stability derivative for three repre- 
sentative airplanes. 
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SYMBOLS 

lateral acceleration, ft/sec^ 
wing span, ft 
trim lift coefficient, 


W 


ip^S 


rolling-moment coefficient. 


Rolling moment 

Jp^Sb 


’hr 


acj 

a** 

2V 


'l o = 


dC, 


'n 


Yawing moment 

^o^Sb 


2V 


_ SC* 

* ^ 
2V 


yawing-moment coefficient, 
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55 * 


c y 






Or = 
r yb 

2V 


lateral-force coefficient, Ia^gral forc e 




D 



°2 


differential operator, d/ds 
Dutch, roll root 

real portion of Dutch roll root 
imaginary portion of Dutch roll root 
root of roll-subsidence motion 


d 3 


root of spiral motion 


fv£\ fa) 
U/l'\P k } 



ratio of nondimens ional rolling velocity to angle of 
sideslip in Dutch roll, roll-subsidence, and spiral 
modes, respectively 





ratio of nondimentional yawing velocity to angle of 
sideslip in Dutch roll, roll-subsidence, and spiral 
modes, respectively 


h 


altitude, ft 


Kn 


constants, functions of measurable flight qu an tities 
(n varies from 1 to 39) 
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nondimens ional radi us of gyration in roll about longitudinal 


stability axis. 


'k- 


: X,0 


2 

COS T) + 


^,0 

b 


sin 2 Ti 


nondlmens ional radi us of gyration in yaw about v ertical 
stability axis, cos ^ + ^ sin ^r\ 

nondimens ional product-of -inertia parameter, 


^Z,0 
b 


E X,0 


sin t| cos T] 


radius of gyration in roll about principal longitudinal 
axis, ft 

radius of gyration in yaw about principal vertical axis, ft 


mass of airplane, slugs 

rolling angular velocity, d//dt, radians /sec 
yawing angular velocity, dijr/dt, radians /sec 
wing area, sq, ft 

nondlmensional time parameter, based on span, Vt/b 

time, sec 

airspeed, ft/sec 

lateral velocity, ft /sec 

weight of airplane, lb 

stability coordinate axes (defined in fig. l) 
angle of sideslip, radians 

inclination of 'principal longitudinal axis of airplane 
with, respect to flight path, positive when principal 
axis is above flight path at nose, deg 

relative density factor, m/pSb 
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P 

4 

♦ 


mass density of air, slugs /cu ft 
angle of roll, radians 
angle of yaw, radians 

DEVELOPMENT OF METHOD 


Equations 

The present method for determining lateral stability derivatives 
from flight measurements depends on arranging the lateral equations of 
motion in such a form that information from each of the three lateral 
modes of motion may he utilized. In this section the original equations 
of motion are given in standard form and then in modified form. 

Assumptions . - The usual assumptions of lateral stability theory are 

made: 

(1) There are only three degrees of freedom: sideslip 3, roll 0, 

and yaw ^ 

(2) After the airplane is disturbed, all aerodynamic controls are 

fixed in their trim position 

( 5 ) The disturbed motions are assumed to be small. 

Mbdlfi cations of standard equations of lateral motion .- The linear, 
second-order, simultaneous differential equations of lateral motion 
referred to stability axes (see fig. l) for the condition of controls 
fixed in level flight are as follows: 

Side force: 

p(2pD - Cy p ) + 0(-C L - JCy p D) 4- t(2pD - Jc^d) = 0 (l) 

Rolling moment: 

P ("%) + ^(^X 2 D 2 - Jc lp D) + ^-2pK xz D 2 - Jc^dJ = 0 


( 2 ) 



6 


NACA TN 40 66 


Yawing moment: 

P(“Cnp) + ^(-2 (iKxZ d2 “ |% D ) + ^K-^D 2 - ^pj = 0 (3) 

Equations (l) to (3) are expressed in a form best suited for deter- 
mination of the resultant motions (3, and ijr) from known system 
constants. The present paper is concerned with the converse problem, 
that of determination of the "system constants (stability derivatives) 
from measured resultant motions. In this procedure, the known quantities 
are the values of period and damping, or time constants, of the various 
modes of motion expressed mathematically as certain values of the 
operator D and ratios of the measured quantities 0 , y, and 3 in 
each mode of motion. If equations (l) to (3) are each divided by 3 
and the terms regrouped, the equations are obtained in a form suitable 
for the present analysis. 

Side force: 

-C Y3 " l^p ^ " 2°*r + + 3 ) - C L | = 0 00 


Rolling moment: 

~ C *3 " i% if " ^ + ^ K xz) = ° (5) 

Yawing moment : 

'% - ^ ^ + 2md(^ *£ - f %z) = 0 (6) 


Basis of the Method 

The present method is based on certain" relations developed in the 
theory of linear differential equations with constant coefficients. These 
relations are proved in most textbooks on differential equations. 

The operator D may be handled as an algebraic quantity in the 
equations of motion. In the usual process of solving these equations, 
the roots of the characteristic equation are obtained (ref. 4 ) . A 
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complex value of D corresponding to one of these roots represents an 
oscillatory mode, the frequency of 'which is related to the imaginary 
part of D and the damping to the real part. A real value of D 
represents an exponential convergence or divergence. In the case of 
the lateral equations of motion, a small real root represents the spiral 
mode, a large negative real root represents the roll-subsidence mode, 
and a complex root represents the Dutch roll mode. 

In each mode of motion, the variables ijr, and p maintain given 

ratios. In the case of nonoscillatory mode, these ratios are real quanti- 
ties. In the case of oscillatory modes these ratios are complex quantities 
the moduli of which give the amplitude ratios between the variables, and 
the angles of which give the phase angles. These ratios may also be 
treated as algebraic quantities in equations ( 4 ), ( 5 ) > and (6). This 
procedure was pointed out in reference 5* 

If the airplane is disturbed in any manner and the controls returned 
to the trim position, the transient motion may be expressed as the sum of 
contributions of the three modes of motion. Because the equations of 
motion are linear, the principle of superposition applies. This principle 
states that, if any transient responses which satisfy the equations of 
motion are added, the resulting response also satisfies the equations of 
motion. Since any response which satisfies the equations of motion con- 
sists of contributions of the three modes, each of these modes must 
separately satisfy the equations of motion. This result is used in the 
present analysis by substituting into equations ( 4 ), (5), and (6) the 
measured values of certain characteristics of each mode. In this way, 
a number of simultaneous equations are obtained from which the stability 
derivatives may be obtained. The details of the procedure are now 
described. 


Method of Solution 


Tabulation of known and •unknown quantities in the equations . - It is 
possible to make ground measurements of all pertinent mass characteristics 
of an airplane. (See ref. 6 .) Furthermore, the air density and the trim 
lift coefficient are easily obtained. Theoretically, all the character- 
istics of each mode of motion could be measured in flight. This informa- 
tion includes the value of the root D and the ratios $/p and ^-/p or 
some derivative of these quantities. These values may be substituted into 
equations ( 4 ), ( 5) > and ( 6 ) to yield three equations for each nonoscil- 
latory mode. The Dutch roll mode yields six equations inasmuch as the 
real and imaginary parts of the equations must equal zero separately. 

Under these conditions, 12 equations could be written to determine the 


nine unknown stability derivatives, namely, Cy , 


P 


C] V 




Cy , 
P 


-v 






and C- 


Because the number of equations available is 
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greater than the number of unknowns , several alternative procedures may 
"be used for determining the unknowns. Some of the equations may be omitted 
from the analysis or some of the quantities measurable in flight describing 
the modes of motion may be eliminated. In practice,, not all the flight 
measurements can be made with equal ease or accuracy. Elimination of the 
measurements which yield the least accurate results would therefore appear 
to be desirable. For the present, one particular set of measurable quanti- 
ties is assumed in order to illustrate the method of solution. Further 
investigation to determine the set of measurable quantities which yields 
the most accurate results would be desirable. The assumption is made that 
all the characteristics of the Dutch roll mode will be used but that, for 
the roll-subsidence and spiral modes, only the roots and not the ratios 
of the variables will be utilized. The following quantities are considered 
as knowns and unknowns: 


Type of measurement 

Known quantities 

Unknown quantities 

Dutch roll 

^r^i'WfOl' 0VP)l 

None 

Roll subsidence 

d 2 

(D9f/j3) 2> (D^/3) 2 

Spiral motion 

D 3 

(d0/p) 3 ,(d t/p) 3 

General 

k xW* k xz^ c l 

None 

Stability derivatives 

None 

VV c v°v 

c n p > ^ l-pf c Y r > c n r > 

c 2 

t r 


Thus, there are 13 unknowns and 12 equations. 


Solution of equations .- Since in the case under consideration the 
number of equations is less than the number of unknowns, some assumption 
is necessary to reduce the number of unknown's. Experience has shown that 
the values of Cy^ and Cy^ may often be neglected in calculating the 

lateral motion. It is therefore assumed in the present analysis that 
these quantities equal zero. A calculation to show the errors introduced 
by this, assumption is given in a later section. 

Substitution of the known quantities into equations (4) to (6) yields 
12 relations, but one of these, the side-force equation for the Dutch roll 
imaginary quantities, is neglected Inasmuch as it involves only known 
quantities if Cyp and Cy-j, are assumed to be equal to zero. The other 
11 equations suffice to determine the 11 remaining unknowns. The method 
of solving the equations is arbitrary, but the procedure outlined in the 
following section has been found to be convenient in practice. 
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For the purpose of the present report, actual flight measurements 
were not used to obtain the values of the "known quantities" listed in 
the preceding table. Instead, numerical values were chosen for the sta- 
bility derivatives, mass parameters, and operating conditions of three 
representative airplanes. These values are listed in table I. These 
selected values were then used to calculate the characteristics of the 
resultant modes of lateral motion. These characteristics, given in 
table II, were then used to supply the "known quantities" listed in the 
preceding table and were used as the starting point in the determination 
of the derivatives by the present method. This procedure assures that a 
consistent set of known quantities are used and, in the subsequent error 
analysis, allows study of the effects of errors in one measurable quantity 
at a time. The practical problems of developing procedures for making 
flight measurements and for working up the desired char ant eristics from 
flight data are not discussed in the present report. 

The procedure for solving the 11 simultaneous equations is now dis- 
cussed. Usually the equations will be solved numerically rather than in 
symbolic form. A numerical example is given in the appendix. In the 
following equations, algebraic combinations of known quantities have been 
represented by (n varies from 1 to 39 ) for the sake of brevity. 

These K values have been worked out in terms of the known quantities, 
but because the resulting expressions are long and because they are not 
needed in a numerical solution, they are not presented here. The fol- 
lowing procedure in terms of algebraic symbols should be used simply as 
an outline of the order to he followed in solving the equations. The 
numerical example given in the appendix shows the actual procedure and 
illustrates how the numerical results may be obtained without specific 
expression for the K values . From the Dutch roll mode, the following 
relations are obtained: 

Equations for imaginary part of Dutch roll mode: 

‘ 2°Ip - 4°lr K l + ^ - 0 <T) 

- + % = 0 ( 8 ) 

Equations for real part of Dutch roll mode: 

+ K 5 = 0 (9) 

~ C h " 2 %^ ' 2 °^*? + ** = 0 < 10 > 

-Cnp “ “ |°Br Kl0 + K 11 = 0 


(H) 
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From these equations, the derivative Cy^ is found and the values of 


, '-'npj w p , aiitt may he determined in terms of Cj r and 


Clp, C np , Cz p , and C 
as shewn in the appendix. 


The solution of the Dutch roll equations is very similar to the 
vector method of reference 3* In the method of reference 3j values of 
two of the derivatives such as and are assumed. The remaining 

derivatives may then he determined by plotting vector diagrams representing 
the equations of motion for the Dutch roll mode. The closure of the vector 
diagrams is the graphical equivalent of setting the real and imaginary 
parts of the Dutch roll equations equal to zero. In the present method, 
hy using the roll-subsidence and spiral as _well as the Dutch roll character 1 
istics, al 1 the important derivatives are evaluated. 


The equations from the roll-subsidence mode may he put in the fol- 
lowing form. In these equations, the K- values are functions only of the 
known quantities, and the previously derived relations have been used to 
express Cj , C^, C^, and in terms of Ci r and C^: 

[¥) z + * Ki ? ■ ° 


(f ) £ ( k 14<\ + k 1 5 ) + (^) 2 + K ie) + K 17°2 r + K l8 - 0 (13) 


(?) 2 (s + Kl9 ) + + K£1 ) + *22% + k 2 }‘° W 


The values of 


(¥), 


and 


a may he eliminated from equations (12), 


(13), and (14) to obtain a relation between and 


c Z r + + K25 = 0 


(15) 


Likewise, the equations from the spiral mode may be put in the following 
form: 
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(^) 3 + *26^)5 + ^7 = 0 (l6) 

(^) 3 ( K 28 c l r + *29) + (|^j 5 ( c I r + k 3o) + K 31 c z r + k 32 = 0 

(Cn p + K33) + (t L ) 3 ( K 3^ P + K 35 ) + K 36 c n p + K37 = 0 (18) 

may be eliminated from equations (l6) , 


fer 


The values of I — I and 

\ P/3 

(17) , and (l 8 ) to obtain another relation between C^ r and 

c Z r + * 38 % + K 39 = 0 (19) 

Equations (15) and (19) may be solved simultaneously to obtain Cnp 
and C-j • Finally, the derivatives C^, C n ^, C-^, and may be 

determined from the equations derived from equations (7) to (ll) . 


DISCUSSION 


Probable Errors in Measurements 

The entire discussion to this point has been baBed on the premise of 
precisely determinable values. Since flight research at best is an 
inexact science, some knowledge of the errors likely to be involved in 
the measurement of the so-called "known" quantities appears to be in 
order. The errors depend, of course, on the accuracy of instrumentation 
employed, the reading accuracy, and the technique of making the flight 
measurements . 

The spiral root is usually the most difficult to measure. A method 
of measuring the characteristics of the spiral mode is given in refer- 
ence 7. The spiral root for a fighter airplane, for which the time to 
double amplit ude was about 30 seconds, was measured to an estimated 
±9-percent error; however, it would be expected that, as the root 
approaches zero, the percent error would increase. Errors in the measure- 
ment of the spiral root are expected to be more of an absolute type of 
error than a percentage error. 

In, order to make a preliminary assessment of the results, the fol- 
lowing ranges of errors, based on flight experience with the type of 
instrumentation employed by the NACA, may be used as a rough guide. 
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Quantity 

Error 

Dutch roll period* percent 

±5 

Dutch roll damping* percent 

±3 

Amplitude ratios* ( Di|r/ 0) and (D0/0) 2 _, percent 

±5 

Phase angle of (Dilr/0)-^ and (1)0/0)^* deg 

±6 

Roll-subsidence root* percent 

±6 

T}* deg 

±1 

Mass and inertia characteristics* percent 

±2 


The mass and inertia characteristics can he measured accurately on 
the ground for a known loading condition* but in flight they are subject 
to uncertainty because of errors in the measurement of fuel consumption. 


In order to determine the effect of errors in the flight measurements * 
calculations were made for three airplanes of low* medium* and high rela- 
tive densities. The stability derivatives* operating conditions* and mass 
parameters which were chosen are typical of the particular class of air- 
plane. These values are listed in table I. These selected values were 
then used to calculate the characteristics of the resultant modes of 
lateral motion. (See table II.) 


Each flight or ground measurable quantity in turn waB varied from 
the correct value and used to recalculate the stability derivatives by 
the method given. The results of these calculations are shown in fig- 
ures 2 to 8. Thus, for example, in figure 2* a 5-percent error in p 
alone (with all the other measurable flight quantities correct) causes 
Cv to be calculated as -0.73 rather than -0.69 and. Ch to be calculated 
*0 T? 

as -0.03 rather than as -0.025 for the airplane with the low relative 
density. 


Errors in Derivatives 


Static stability derivatives . - Examination of figures 2 to 8 shows 


that the static stability derivatives C- 


Y n ' 


p' V 


and C n ^ are* in general* 


not unduly sensitive to errors in the measured data. The percent error in 
these derivatives is usually about the same or less than that in a given 
measured quantity. The errorB usually increase somewhat as the value of 
p increases. One exception is the relatively high sensitivity of Cy 

to the phase angle of Dijf/0 for the Dutch roll mode (fig. 6) . 


■0 


A surprising result is the insensitivity of C n to the frequency of 


0 


the Dutch roll (See fig. 7.) The frequency is usually considered 

to be proportional to the square root of C n for constant values of the 
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other derivatives. In the present example, the large simultaneous varia- 
tions of the other derivatives, especially C n ^, probably account for the 

small variation of Cn R . It is not known whether this result would be 
found in all cases. * 

Damping derivatives .- The damping derivatives Ci and C n become 

progressively more difficult to determine as the value of p increases. 
This result is related to the fact that the damping derivatives have leBB 
effect on the motion as p increases. Accurate determination of these 
derivatives is therefore less important at higher values of p. For the 
range of p considered, reasonably accurate measurements of were 

obtained. The value of Cj^ can be determined only for the airplanes 

with low and medium relative density. 

Cross derivatives . - The cross derivatives C? and C- are evi- 

6 r °p 

dently the most difficult to obtain accurately. The value of C, is 

t r 

well-defined for the airplane with the low relative density but rather 
poorly defined for the airplanes with medium and high relative density. 
The value of is not determined with reasonable accuracy for any of 

the cases. It is greatly affected by errors in Dutch roll period and 
in the phase angles of D^/p and DtJt/{ 3 . 


Effect of Side-Force Rate Derivatives 


A check was made to determine the effects of variations In the side- 
force rate derivatives Oy and Cy^ which were assumed to equal zero 

In the previous analysis. 


r ~P 

In practice, the values of 


may 


Cy^ and Cy^ 

lie in the range from about 0.3 to -O.3. The values of all the other 

derivatives were calculated with these values for Cy anrl Cv . Tha 

r x p 

results of these calculations are shown in figure 9. The derivative Cy 

x r 

has very little or no effect on any of the other derivatives. The deriva- 
tive Cy^ has a small effect on Cy^ and C^ for all three airplanes 

and a small or no effect for the other derivatives. In view of these 
results it appears that a slight improvement in accuracy could be gained 


by -using an estimated value for 
assuming it to be zero. 


Cy In the analysis rather than by 


Effect of an Error In p 

The separate errors assumed In the quantities D 0 /p and Dijr/p may 
be considered as errors in D$ and Dtjr. If it is assumed, however, that 
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an error exists in the magnitude or phase angle of (3, "both D$/p and 
D^/p will "be affected in the same way. Additional calculations have been 
made, therefore, to determine the effects of the combined errors in these 
quantities which would result from an error in p and the results of 
these calculations are given in figure 10. The derivative Cy. is highly 


sensitive to the phase angle 
tude of p. The derivatives 


P 

p hut shows little sensitivity to the magni- 


'V 


J v 


and CL become progressively 


more sensitive as the value of u, increases whereas C 7 , C 7 , and 

V V 

for all three airplanes show very small or no sensitivity to p. 



Improvements in the Method 

Measurement of Cy^.- A deficiency in the method described is the 
sensitivity of Cy^ to an error in the phase angle Ity/p. (See fig. 6.) 

The source of this error may be seen by examining the side-force equation 

for the Dutch roll mode. (See eq. (4).) Since Cy_ is the only unknown 

P 

involved, this equation alone is required for its determination. 


-% + ^ + - °4l,r - ° 


( 20 ) 


Since Dty/p is almost a pure imaginary quantity (see table II) , changes 
in the phase angle affect the real part of this quantity directly. This 
value is then multiplied by the large factor p which introduces errors 

in c Y p - 


An alternate method of determining Cy , pointed out in reference 5, 

P 


is to measure lateral acceleration a_, by means of an accelerometer 

y /_ \ /.\ 


C L 


which 


The value gives the sum of the terms 2uD + 2u.[2£) - 

PqS * * U / 

according to equation (20) equals Cy • In addition, the use of this 

P 

procedure may allow determination of Dijr/p to a higher degree of accuracy 
than would be possible by measurements of Dijr and p directly. 


Alternate sets of flight measurements .- The set of flight measurements 
discussed in the present report is not necessarily the best for any particu- 
lar airplane. Since determination of C n ^ by the present method requires 

extreme accuracy in the measurements, it is likely that the use of addi- 
tional data, such as the value of Dy/p for the roll-subsidence mode, 
would allow more accurate determination of C n . Preliminary calculations 
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indicate that greater accuracy in the measurement of Cnp can he obtained 
by this method. 

Another possibility which may prove desirable in some cases is elimina- 
tion of measurements of the spiral mode and substitution of more complete 
data on the roll-subsidence mode. Accurate determination of the spiral 
characteristics is difficult in some airplanes because of the changes in 
lateral trim caused by motion of fuel in wing tanks or because of the 
tendency of wings to take a slight twist due to structural hysteresis 
following application of lateral control. 


CONCLUDING REMARKS 


A method is presented for determining lateral stability derivatives 
from flight measurements . This method utilizes data from each of the 
three modes of lateral motion and allows determination of all the important 
stability derivatives. An error analysis is made to show the effects of 
errors in the measured quantities on the accuracy of each stability deriva- 
tive for three representative airplanes. 


The static stability derivatives Cy , 

3 


C 7 , and 

V 


n 


3 


may be deter- 


mined with good accuracy by the proposed method. The damping deriva- 
tives C^p and C nr become progressively more difficult to determine 

as the value of the airplane density factor p increases. The values of 
the cross-rate derivatives C^ and C^ are the most difficult to obtain 

accurately. The value of C^, is well-defined for low values of p but 

is poorly defined at medium or high values of p. The value of C n ^ 

is not determined with sufficient accuracy for ary of the cases studied. 


The proposed method may be modified to utilize different sets of 
flight measurements without changing the basic procedure. Improvements 
In the accuracy of determination of some of the derivatives might result 
from the use of different sets of flight measurements. 


Langley Aeronautical Laboratory, 

National Advisory Comnittee for Aeronautics, 
Langley Field, Va., May J, 1957- 
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APPENDIX 

NUMERICAL EXAMPLE 


The data given for the fighter airplane in tables I and II are used 
in the example. Values are substituted into equations (it -) , ( 5 ), and ( 6 ) 
of the text. 


Imaginary Part of Dutch Roll Mode 


The side-force equation disappears because the two unknowns 

and Cy are assumed to be equal to zero. The rolling-moment equation 
i r 
is: 




- Jc Zp (0.103i) - Jc Zr (- 0 . 302 i) + 2(13) (-0.035^ + 0.3041) [(- 0.211 + 
0 . 1031 ) ( 0 . 0171 ) - ( 0.0100 - 0 . 3021 ) ( 0 )] = 0 


In the products (O.O 354 + 0. 304i) (-0.211 +0. 1031) and (-0.0354 + 
0. 304i) (0.0100 - 0 . 302 i), only the Imaginary parts are vised. 


- 0 . 0515 iC 


h 


+ 0 . 15 UC Zr + 0.03021 = 0 


C ip = 2.90 C^ - 0.585 


(Al) 



p 
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The yawing-moment equation is: 



- Jc. (0.1031) - Jcn^-o.3021) + 2(13) (-0.0354- + 0.304-i) [(0.0100 - 

0.302i) ( 0 . 04 - 92 ) - (-0.211 + 0.103i)(0)] = 0 

Again, in the products (- 0 . 0354 - + 0 . 304 -i) ( 0.0100 - 0 . 3021 ) and 
(- 0 . 0354 - + 0 . 304 i) (-0.211 + 0.1031) , only the imaginary parts are used. 

-0.051510^ + 0 . 15110^ + 0.01751 = 0 

Cn r = 0 . 340 Cnp - 0.116 (A 2 ) 


Real Part of Dutch Roll Mode 
The side-force equation is: 



If the following equalities are used, the substitution in the side-force 
equation may he followed with ease. 
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Substituting the values in the side-force equation gives: 


Cy + 0 + 0 + 2(13) (-0.0354) + 2(13) (0.0100) 

0 


- °- 0T {% 


- 0.211 


0354 


+ 0 . 103 i \ = 

+ 0 . 304i J 


= 0 


In the quotient of 


-0.211 + 0.1031 
-0.0354 + 0 . 304i 


) 


only the real part is used; thus, 


-Cy^ - O .69 = 0 


Cy = -0.69 


(A3) 


The rolling-moment equation is: 

■% ■ 5%(ir) ljr ■ + ^ 


T 


2 I £* 
P 


J K X " l 3T K XZ 


= 0 


-C^ " |%(- 0 - 211 ) " |c lr (0.0100) + 2(13) (-0.0354 + 0.3041) [(-0.211 
0.103i) (0.0171) - (0.0100 - 0.3021) (0)J = 0 

In the products (-0.0354 + 0. 304i) (-0.211 + 0.103i) and (-O.O 354 + 
0*304i) (0.0100 - 0 . 302 i), only the real parts are used; thus, 

-C^ + o.io6c z - 0.0050^ - 0.0105 = 0 

Substituting the value of from equation (Al) gives: 

-Cj + 0. 106^2. 90C Z ^ - 0.585) - 0.0050^ - 0.0105 = 0 


Cj = 0 . 300c ^ - 0.0720 


(A4) 
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The yawing-moment equation is: 

'S ' 4sWj,r ' + ^ 


'Sat) k/ . fall k. 

VP/I z \p 


■xz 


= 0 


-CL. - icL (-0.211) - JCL, (0.0100) + 2 (13) (-0.0354 + 0-304i) (0.0100 - 

u p 2 2 L 

0.3021) (0.0^92) - (-0.211 + 0.103i) (0)J =0 

In the products (-0.0354 + 0. 30lj-i )( 0.0100 - 0.3021) and (-0.0354 + 
0.304i) (-0.211 + 0.1031), only the real parts are used; thus, 

-c no + 0.1060^ - 0.0050^ + 0.116 = 0 


°r 


Substituting the value of C n from, equation (A2) gives: 

- c n p + 0.1060^ - 0.005(0.340C np - 0.116) + 0.116 = 0 


c nm = 0.1050^ + 0.117 


°P 


♦ 


Roll Subsidence 
The side-force equation is: 


— Cy - ^Cy feA - ^Cy + 2pD 2 

^P 2 Y p\p ^ 2^r\^p j 2 ^ 2 


1 + 'f 


^ C T = 0 


(A5) 


0.69 + 0 + 0 + 2(13) (-0.499) + 2(13) 



2 



= 0 



(A 6) 



20 


NACA TN 4066 


The rolling-moment equation is: 


-% - i°h.(?) P ■ + ^ (tL ** 2 ‘ (?„ 


K- 


XZ 


"% " 4%(l^) + 2(l3)(-0.499) (0.0171) (0) 


Substitutions from equations (Al ) , (A4 ) , and (A 6) give: 


-(0.5000^ - 0.0720) - ^2.90Cj r - O.585) 


**f) - 86.5 


|o, fst) + 2(13) (0.499) 

2 & r s p y 2 


185 (a) . 86.5 


(0.0171) = 0 


(l 2.8 - 265 C, V^) - 5.95 + I 25 C 7 «= 0 

V r /\P/2 r 


-5.95 + 125C, 


m = 

Py 2 -12.8 + 2650, 


(A7) 


The yawing-moment equation is : 


-c n - iQn - icvf^ + 211D2 

n 0 ^VP/o 


^ Kz 2 - (& 


KXZ 


"S ‘ ' IS^) + 2(13) (-O-^) (^) (0.0492) - ( 0 ) 


= 0 
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Substituting the values obtained from equations (A2 ) } (A5 ) , a nd (a 6) 
gives : 


-(° 


.1030^ + 0.117) - ^ 





^ 0 . 3400 ^ - 0.116 ^ + 


(2) ( 13 ) (-0.499 ) {^j (0.0492) = 0 




+ 43. 


3 % “ 


0.117 = 


0 


■ p*\ -0.117 4-43.30^ 

,P/ 2 ~ 0 . 58 + 91 . 7 Cnp 


(A8) 


Substitute the value of 


(lty/f3) 2 from equation (A7) into equation (A8) . 


-5-95 + 125 C -0.117 + 43-50^ 

-12.8 + 2650^ 0.58 + 91.70x1^ 


103.50, + 8. 

h r 


65 %- 


4.95 = 0 


(A9) 


Spiral Motion 


The side-force equation is: 



' 2%(f \ - + ^ _ X + (|) 5 



= 0 


0.69 + 0 + 0 + 2(13) (-0.0000725) + 2 (l 3 )f^l - ( Q ‘° 71 W) = 0 

V p /3 V-o. 0000725 /\P j 3 

2 6^ + 980^ + 0.688 = 0 
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= -37. 7^) - 0.0265 (A10) 

\fi) 3 W 3 


The rolling-moment equation is: 


-C 


Z 3 



" 


Di|f 

3 





Substituting the values obtained from equations (Al ) , (A4 ) , and (A10) 
gives : 



^i7.oc 2r + 0.292 - 0.3000^ 


+ 0.0720 = 0 



0.0720 - 0.3000^ 
-0.292 - 17.00^ 


3 


(All) 
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The yawing-moment equation is: 


-s - 1 ■ fe ^ 2 - (f) 


E 


■xz 


= 0 


"S - 2\(f) 3 - &r($ 3 + 2(13) (-0.0000725) (^(0.492) - (^(0)| 
Substituting the values from equations (A2 ) , (A5 ) , and (A10) gives: 

-(0.103^ + 0.117) - “ |(°- 340 C np - 0.116) 37 . 7 ^) - 




■ 

0.0265 

- 0.00188 

- 37 • Tl^) - 0.0265 

- 


VP / 3 


(0.0492) = 0 



2 . 20 ) . 0.099^ - 


0.119 = 0 



-0.119 


2.20 


' °-°99Cnp 
5.400^ 


(A12) 


Substitute the value of (D^/p)* from equation (All) into equation (A12) 
to obtain ^ 

0.0720 - O.3OOC7 -0.119 - 0 . 099 C,, 

r _ iip 

-0.292 - 17.0C , 2.20 - 5.40C,, 


-2.68C, 

4 r 


0.420C„ + 0.124 = 0 


(A 13 ) 


= -0.l57C np + 0.046 
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Substitute the value of C 7 from equation (A13) into equation (A9) 

103.5^0.1570^ + 0.01+6^ + 8.630^ - 4.95 = 0 
-7.670^ - 0.19 = 0 
= -0.025 

Substitute the value of into equation (AI 3 } to obtain 

Gi r « -0.157( -0.025) + 0.046 
c, = 0.050 

t r 

Substitute the value of Cj into equation (Al) to obtain 

c ^ = 2.90(0.05) - 0.585 

C 7 = -0.44 

Substitute the value of into equation (A2) to obtain 

Cj^ = 0.34(-0.025) - 0.116 


= -0.125 


Substitute the value of Cj r into equation (A4) to obtain 
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Substitute the value of 



into equation (A5) to obtain 
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TABLE I 

CHARACTERISTICS OF REPRESENTATIVE AIRPLANES USED IN CALCXJLATIONS 


Char act er is t I c 

Fighter 

Medium Bomber 

High-altitude fighter 

n 

13.0 

31.83 

182 

% 2 

0.0171 

0.0311 

0.01557 

tP 

ro 

0.0492 

0.072 

0.156 

K xz 

0 

0 

0.0020 

c L 

0.071 

0.443 

0.49 

h, ft 

0 

35,000 

50,000 

V, ft/sec .... 

700 

700 

776 

h, ft 

41.6 

116 

25 

Cy 

-0.69 

-0.61 

-O.58 

Cy 

Z P 

0 

0 

0 

Cy 

i r 

0 

0 

0 

*p 

-0.0573 

-0.14 

-0.18 

% 

-0.44 

-0.44 

-0.33 

C lr 

O.05 

0.149 

0.22 

% 

0.H5 

0.12 

0.25 

% 

-0.025 

0.0276 

-0.049 

Sr 

-0.125 

-0.156 

-0.68 







ro 

00 


1ABIX II 


GAICnUUBD CHABflraffllSffilCS OF THE EESOI/EMH? M 3 DE 8 HK 1HE REFKSSHfEffrrVS AIHPLAHffl 


j^All roots are given in ncradimnsional farm.^] 


Airplane 

li 

Roots of characteristic equation 


D*/P 

Dutch roll 

Boll 

Bubflldence 

Spiral 

notion. 

Dutch roll 

Roll 

gnbsldance 

Spiral 

motion 

Dutch roll 

Roll 

subsidence 

Spiral 

■otton 

7 igtrt«r 

15.0 

- 0.0354 ± 0.30591 

-0.4993 

-O.OOOO725 

-0.2U3 + 0,10281 

24.77 

- 0.04947 

0.01005 - 0.30221 

0.3375 

1.84 

Madluu 

boober 

31.83 

- 0.00447 ± 0 .l£ 79 i 

-0.1284 

- 0.000419 

-0.215 + 0.28281 

4.36 

- 0.095 

0.00684 - 0.1591 

-0.1177 

1.56 

High- 

altltula 

filter 

182 

0.00258 t O.O665I 

- 0.0410 

-0.000770 

-O.I97 + 0.37451 

2.75 

- 0.49 

O.OO325 - 0.06221 

-0.0508 

0.836 
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Amplitude error, percent 


Phase angle error, deg 


Figure 6.- Variation of the stability derivatives produced by errors in 
the amplitude and phase angle of the Dutch roll ratio Dijr/p. 






Error ia Dj, percent 


Error in D3, per cant 


Figure 8.- 
the roll 


Variation of the stability derivatives produced by errors in 
•subsidence and spiral roots of the characteristic equation. 








